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Abstract 

Propagation of P and SV waves in an elastic solid containing randomly dis- 
tributed inclusions in a half-space is investigated. The approach is based on a mul- 
tiple scattering analysis similar to the one proposed by Fikioris and Waterman for 
scalar waves. The characteristic equation, the solution of which yields the effective 
wave numbers of coherent elastic waves, is obtained in an explicit form without the 
use of any renormalisation methods. Two approximations are considered. First, for- 
mulae are derived for the effective wave numbers in a dilute random distribution of 
identical scatterers. These equations generalize the formula obtained by Linton and 
Martin for scalar coherent waves. Second, the high frequency approximation is com- 
pared with the Waterman and Truell approach derived here for elastic waves. The 
Fikioris and Waterman approach, in contrast with Waterman and Truell's method, 
shows that P and SV waves are coupled even at relatively low concentration of scat- 
terers. Simple expressions for the reflected coefficients of P and SV waves incident 
on the interface of the half space containing randomly distributed inclusions are 
also derived. These expressions depend on frequency, concentration of scatterers, 
and the two effective wave numbers of the coherent waves propagating in the elastic 
multiple scattering medium. 
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1 Introduction 

We consider the problem of elastic wave propagation in heterogeneous solids containing 
distributions of inhomogeneities. The results have application in geophysical exploration 
and ultrasonic evaluation of composite materials or biological tissues. Typically, the 
inhomogeneities can be hard grains, inclusions, micro-cracks, fibers, pores, or contrast 
agents, for instance. The present work is a contribution in this direction. 

The main difficulty is the quantitative description of interactions between many scat- 
terers distributed throughout the heterogeneous solid. In the case of random sets of 
scatterers this problem cannot be solved exactly, except in very particular cases [25] , and 
in general only approximate solutions are available. These are all based on hypotheses 
that reduce the many scatterer problem to a problem for one scatterer or one cell. In this 
way we may focus attention on the coherent wave propagation, which is the statistical 
average of the dynamics corresponding to all possible configurations of the scatterers. It 
is well-known that the coherent motion makes each heterogeneous medium appear as a 
dissipative homogeneous material, and propagation is governed by a complex effective 
wave number that is frequency dependent. The real part of the wave number is related 
to the velocity, and the imaginary part represents the attenuation. 

Various methods for the analysis of the multiple scattering of elastic waves have been 
considered in the literature. These include: perturbation methods [HJ [27]; self consis- 
tent methods, which are analogous to the coherent potential approximation [9j [15]; and 
effective medium methods based on self consistent schemes j26j [5j [13]. The method of 
asymptotic homogenization [21] was recently considered in the context of fibre-reinforced 
media where the fibres are arranged on a periodic lattice. All of these approaches have 
their own advantages and disadvantages. 

The method developed here starts from an explicit multiple scattering formulation in 
which the field scattered from any particular scatterer is expressed as a multipole (far- 
field) expansion. In the case of scalar waves, this is a classical topic with a large literature. 
The method was initiated by Foldy for studying the multiple scattering of isotropic point 
scatterers [11]. The Quasi Crystalline Approximation (QCA) was subsequently introduced 
by Lax in order to take into account double scattering [16]. With the use of the QCA, 
Waterman and Truell derived an expression for the effective wave number of coherent 
waves that depends on the far-field scattering properties of a single scatterer, in both 
the forward and backward directions [33]. Waterman and Truell's approach, like that of 
Twersky [29] which is based on the Foldy's approximation, is known to be valid only for 
dilute distributions of scatterers. Not long after, Fikioris and Waterman introduced the 
"hole correction" that allows consideration of higher concentrations of anisotropic finite- 
size scatterers [10] . Anisotropy refers here to scattering properties that depend on angular 
orientation, as illustrated for example in a polar radiation diagram. The present approach 
is based on multipole expansions, and can therefore handle scatterers with non-circular 
or non spherical shapes. The T-matrix method provides a convenient tool to calculate 
multipole expansions [3"T] . 

Numerous investigators have attempted to improve on the original theory [HI [331 02] 
or extend its range of applicability. Recently, Aristegui and Angel [3] showed that Foldy's 
method also yields predictions of the reflection and transmission of scalar waves by a 
random distribution of point or line scatterers contained within a slab. Similarly, an 
extension of Fikioris and Waterman's formalism was developed by Le Bas et al. [TTJ in 
order to describe both the reflection and transmission from a slab-fluid region in which 
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elastic cylindrical scatterers are randomly spaced. They show in particular that two 
different coherent waves can propagate in the fluid slab. Starting from Fikioris and 
Waterman's approach, Linton and Martin [18] obtained the two-dimensional counterpart 
of Lloyd and Berry's effective wave number [201 [T9]. taking into account the hole correction 
and the boundary effect. In the case of immersed fluid cylinders with a sound-speed that 
is close to that of the surrounding fluid, Martin and Maurel [22] show that Linton and 
Martin's formula [TS] can be derived from Lippman-Schwinger's equation without use of 
the QCA, thus providing a link between distinct formalisms based on Green functions and 
multipole expansions. Shear-horizontal elastic waves (SH) are also scalar-type waves, and 
the effective dynamic properties SH waves in composites made of elastic cylindrical fibers 
randomly distributed in another elastic solid have been calculated by Aguiar & Angel [2J 
and Aristegui & Angel [I]. They show that the effective mass density and the effective 
shear stiffness are complex valued and frequency dependent. 

In comparison with the numerous studies of the scalar situation, multiple scattering 
of elastic waves involving both compressional (P) and shear waves (SV) has received 
relatively little attention. Varadan et al. [30] arid Yang & Mai [31] have considered 
this problem using multipole expansions. Their analysis is mainly focused on the low- 
frequency limit (Rayleigh limit) which predicts dynamic effective mechanical properties 
of particulate composites that are in agreement with Hashin and Rosen's bounds [T2] . 
The Generalized Self Consistent Method (GSCM) [31] is derived by using a self consistent 
scheme applied to Waterman and Truell's formula [23] ■ However, this formula, which is 
valid for scalar waves, is applied to P and SV waves separately [31], and mode conversions 
between P and SV waves are therefore neglected. In contrast to the GSCM, the theory 
developed in [30] takes mode conversions into account, but the equations that involve P 
and SV waves are uncoupled by invoking additional hypotheses above and beyond those 
for the QCA used in most of the papers previously cited. This is not necessary and 
it can be avoided by ensuring that the equations faithfully and accurately describe the 
coupling between the P and SV waves, as we do in this paper. More precisely, we derive 
equations for P and SV waves in the spirit of the paper by Fikioris and Waterman [TO] , 
and effective wave numbers are then obtained in the limit as the scatterer size tends to 
zero. From this point of view, our results can be considered as a generalization of the 
work of Linton & Martin [18] for scalar waves to the elastic case. It is shown here that 
the effective wave numbers, which are predominantly P or SV waves at low concentration 
of scatterers, depend upon the four far-field scattering functions, including those related 
to mode conversions. An advantage of our method is that it is valid not only at low 
frequency but also at higher frequencies. For this reason we also consider the short-wave 
limit in some detail. Also, we assume that the scatterers are located in a semi infinite 
region, a configuration that has received little attention. This allows us to determine the 
four reflection coefficients related to the solid-solid interface, both in general and in the 
small concentration and long wavelength limit. 

The outline of the paper is as follows. Section |2J contains a derivation of the coupled 
equations for elastic waves. This is based on the multiple scattering theory of Fikioris 
& Waterman [TU] for acoustic waves. In Section [3] we reformulate the infinite system of 
equations to derive a single equation for the wavenumber. We present in Section H] a 
systematic solution of the coupled equations using asymptotic expansions in the concen- 
tration. This generalizes the theory of Linton & Martin [THl [TU] . The high frequency 
limit is also derived and compared to the Waterman & Truell approach. Section is 
devoted to the calculation of reflection coefficients at the interface of the half-space en- 
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closing the random distribution of scatterers. Asymptotic approximations are derived for 
the reflection coefficients that are valid in the low concentration limit. In Section |6] we 
derive explicit formulae for low concentration expansions of the effective wavenumbers, 
formulae which generalize that of Linton & Martin for the scaler case. We begin with the 
multiple scattering equations for elasticity. 

2 Fikioris & Waterman theory for elastic media 

The multiple scattering theory of Fikioris & Waterman [10] was derived for acoustic media. 
Here we develop the theory to consider multiple scattering from identical scatterers in 
elastic materials. 

2.1 Multiple scattering equations 

Suppose that time-harmonic P or SV waves are propagating perpendicular to N parallel 
cylinders located in an elastic solid and that and kx are the wave numbers of the P 
and SV waves. We assume the Helmholtz decomposition of the displacement in the form 

u = Wip L + V x (i) T e z ) (1) 

where e z is the unit vector parallel to the cylinders and ip L and i/j T represent potentials for 
the longitudinal (P) and transverse (SV) components of the waves. Under the influence 
of the incident waves ipi nc (r) and ?pf nc (r), both L and T scattered waves ^s^^fc) an d 
IP's {f, fk) are generated by the kth scatterer, so that 




Here, the first vector argument f specifies the field point of evaluation, while ff. is the 
location of the kth scatterer. The scatters are assumed to be identical in composition and 
orientation and the properties of a single scatterer are assumed to be known, so that a 
rule is available that relates the scattered waves ^(r, r k ) and the exciting fields ip%(f, fk) 
acting on the kth scatterer (a = L,T). This rule defines a linear scattering operator 

fips(r,r k )\ = 

yP T s if,f k )) ' 

If the scatterers are similar but are not all oriented in the same manner then the orientation 
of scatterers must be assumed to be random. The exciting field acting on the kth. scatterer 
is the sum of the incident waves and the scattered waves from all scatterers other than 
the kth. It follows that 




Equations (j2J) to (j3J) are the multiple scattering equations that generalize those obtained 
for the scalar case (cf. eqs. (2.9-10) in [SB]). 
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In order to derive the equations governing the coherent motion, we use the method 
initially developed by Foldy pj] and Lax [16] to average over all configurations of cylinders. 
This method is a very well documented [211 EH], an d it includes as a special case the quasi- 
crystalline approximation (QCA). Performing the configurational average transforms eqs. 
(H) into 



<^(r,n)> 
(V>I(r,ri)) 



dfj n(fj, fi) 



T 



LLi 



T 



TLi 



(V4(r,r})> 



(ipl{f,fj)) 



(5) 



In these equations, f\ is the location of one of the cylinders, {^%{f,fj')) 



a 



L, T) are 



the average coherent fields acting on the jth scatterer, n{f,fj) the conditional number 
density of scatterers at f if a scatterer is known to be at r}, and the integral is taken over 
the whole surface accessible to scatterers. For a uniform and random array of identical 
cylinders of constant density n and radius a, the "hole correction" [10J requires 



n(f, fj) 



n for \f — fj\ > b, 
otherwise, 



(6) 



with b > 2a. Generally, b represents the distance of closest approach between centers 
of adjacent cylinders. Equations ()5]) are the multiple scattering integral equations that 
generalize the integral equation for acoustic media (cf. eq. (2.1) in [10]). in the same 
way, the average coherent fields are obtained by performing a configurational average on 
eqs. 02]), yielding 



<^(r1> 



lb] 



t rie \ 



+ n Q / d r 



T LL {f 3 ) T TL (r- 



(^(f,r})> 

(^E^fj)) 



(7) 



It should be noted that eqs. ([7]) are exact if the observation point r is not close to the 
surface of the volume enclosing the cylinders. Field points at the boundary of the volume 
require a special treatment, see [18] for details in the acoustic case. 



2.2 Modal equations 

Our goal in the remainder of this Section is to obtain the modal equations, which will be 
solved in Sections [H] and H] to give the effective wave numbers. In the following, cylinders 
are assumed to be randomly distributed in the half-space defined by S + = {x > 0}. The 
incident wave propagates in the direction of the x— axis, normal to the interface x = 
(r = (x,y)): 

^ nc (f) = A a e^, «G{L,T}, (8) 

with Al,t = 1 or depending on the type of the incident wave considered (P or SV). The 
incident and the effective potentials satisfy the Helmholtz equation 

V 2 f + fc> a = 0, ae{L,T}. (9) 

In the same way as ijij 1 \f,fj), the effective potentials (^e' T (^ ^j')) satisfy the Helmholtz 
equation and are regular functions at the point fj, they can therefore be expressed 

(Mr, fj)) = £ Al{fj) J n (k aPj ) e me ^\ a G {L, T}, (10) 
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with pj = f — fj, 9{pj) = arg(pj) and pj = \pj\. As usual with the T-matrix approach 
[31] . the transition operators are defined by 

T a/3 (r"}) J n (k aPj ) e me ^ = TfH^(k m ) e md ^\ a, (3 E {L, T}, (11) 

and the corresponding far-field scattering amplitudes of the different interactions are given 
by 

T a P(6)e ikaX = ^i n TfH^\k p r) e md 



e^-^/^W, a,/3G{L,T}, (12) 



with f — (r cos 9, r sin 9). The far-field scattering functions f a ^(9) are therefore Fourier 
series with coefficients equal to the modal scattering amplitudes T"' 3 , i.e., 

r?(9) = Y,Tfe in \ a, /3e{L,T}. (13) 

n 

Modal coefficients T^ 13 can be calculated numerically j3TJ [32]. For circular cylinders, 
they are the components of the T-matrix and satisfy the symmetry relation = . 
For non circular cylinders, they are expressed in terms of the T-matrix components that 
depend on the orientation of the scatterer, so that in general ^ . 
Substituting eqs. (fTUl) and (fTTj) into eqs. ^ gives 

J2l A n(ri) -i n A L e tk ^}J n {k LPl )e md ^ 

n 

Js+ 

n 

+ J2 T n L I df j n(f j ,f 1 )AZ(fj)Hi 1 \k LPj )e i ^\ (14a) 

(L^T). (14b) 

The symbol (L-B-T) means the same equation(s) as the previous one(s) but with the L 
and T indices permuted. In order to write series of (|T4|) as a function of the coordinates 
centered on f x = (xi,yi), we use the change of variables 9(pj) = 9(—pj) ± ix with the 
addition theorem pQ 

H^(k a p 3 )e me ^ = ^(-^-^^-^^^^^^^^^(^pOe^^), (15) 

in 

where fji = fj — fi, 9{fj\) = arg(r jl ) and Tj\ = \ fji\. The condition \px\ < \ fji\ must hold 
in order that the waves scattered from fj are expressed with the use of eq. ( I15p in the 
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vicinity of f±. It then follows that 

ifj nfo, ^^foO^'WiO e ip6 ^ +i n A L e ife ^, (16a) 

(16b) 



_->\p T TL 
1 ) ± n+p 



s+ 



(LhT). 



We seek coherent waves that propagate in the equivalent homogeneous medium, which 
consists in assuming that solutions of eqs. (fTBT) may be written in the form 



r a l 




( e *- r * 1 


= i n 








A T 










Bl_ 





(17) 



Here £ = (£, 0) and £ = (£ , 0) are the effective wave vectors of coherent waves that 
propagate in the direction of the x— axis, and the coefficients A^ and B% are at this 
stage unknown. Here, two coherent waves with £ and £ as wave numbers are assumed 
to propagate, which is a natural hypothesis for scarce concentrations of scatterers. In 
such situations the homogeneous medium looks like an elastic medium in which the two 
waves that propagate are predominantly P or SV waves (see §4. 1 j) . In the following, the 
concentration is supposed to be low enough that only two coherent waves propagate. It 
is worthwhile to note that a search for solutions of eqs. ([TBI in the following form 



(A L n (r 3 ; 



A L n e* 



A L 

Al^ 



(18) 



leads to an ill-conditioned problem (the number of unknowns does not equal the number 
of equations available). Solutions of the form (JTSJ) were considered in [30], although in 
the context of uncoupled integral equations. 
Inserting the ansatz (JT7J) into eqs. ([TBI) yields 



A L n +B% e^'- rl = A L e ikLXl + % 



with 



T n + p A n+p I p (£) + T n+p B n+p I p (£ 

, t TLaT jL( c\ _j_ rpTL tz>T tL(c;'\ 
~^ 1 n+p I± n+p 1 p \S>) 1" 1 n J rV D n j rp l p \Z, ) 

(L++T). 



df j n(f J ,f 1 )H p :i \k a r 3l )e ip9 ^e i ^, ( E {L,T}. 



W) = 

Taking the hole correction ([6]) into account, the integral ( 1201) may be written 
lp(0 = n e^ [ dr^MKr^e 1 ^, ( E {£,£'},«€ {L,T}, 



(19a) 
(19b) 

(20) 
(21) 
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where Sj, is the entirety of S + with the exclusion of the hole \ fji\ < b centered at r*i, and 
<pp are the cylindrical wave functions with respect to the same origin, with 

<f> p (k arjl ) = H^(k arjl )e ip9 ^\ a G {L, T}. (22) 

The integral in eq. f[2"TT) has been previously calculated in [T71 [£B] , with the result 

ff(0 = ^[#^^(C)e^ + e^], Ce {£,£'}, «G{L,T}, (23) 

where 

^(C) = C^(C6)^ 1) (M-W(C6)^ 1) '(A; a 6), Ce {£,£'}, ae{£,T}. (24) 
As a consequence eqs. ( IT9|) become 

A L e «xi +B L gi^xi gifc £ *i 



^ c t2 _ £.2 / / L n+p^rt+p ~ 1 n+p^n+pi 1 y p K'i) 
S L p 

> L p 

l pifc^! 2m \^[T LL 4 L 4- T TL 4 T 1 



(L-B-T). (25b) 
2.2.1 The generalized Fikioris & Waterman equations 

Equations (j25j) are satisfied whatever the value of X\ if the coefficients of e^ Zl , e^ x i ; e lfciXl 

and e tkTXl are all zero. Equating the coefficients of e*^ 1 and e 1 ^ Xl to zero gives what is 
known as the Lorentz-Lorenz law. We obtain 

A " ~ E< T "4 + T p LA p) N n- P (0 = 0, (L^T), (26a) 

S L p 

B L n - £ (3^^ + = 0, (L*»T). (26b) 



Equations (1261) are the desired modal equations that generalize those of Fikioris and 
Waterman. They provide two identical homogeneous systems of linear algebraic equations 
which involve either the unknowns {A^A^} with £ or {B^B^} with £ . The existence 
of nontrivial solutions of the homogeneous system determines the effective wave numbers 
£ and £ . This is the goal of the next Section. 
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2.2.2 The extinction theorem for elastic waves 

Equating the coefficients of e lkLXl and e tkTXl in ( 125|) to zero corresponds to what is known 
as the extinction theorem. In this case it becomes 

Al + *J^k) £ [T ' LLA * + T " A ' ] + dr^i S [r " s * + r " s?] = - 

(27a) 

(L-<->T). (27b) 

The extinction theorem is useful for calculating reflection coefficients for the waves re- 
flected at the interface x = 0, see Section 

2.2.3 The Waterman & Truell method 

Waterman & Truell's approach is an alternative method relevant at very low concentration 
of cylinders [3] and simpler than the more general approach of Fikioris & Waterman 
represented by the system (|26|) . The Waterman & Truell approximation assumes a pair 
correlation function with the following property [33| [18] 



/ — * — * \ J™o for |x - Xi| > 77, 
w(r,ri) = <_ . (28) 

I (J otherwise, 

for 77 — > with r*= (x, ?/) and fi = (xi,?/i). In this limit eqs. (120]) are improper integrals 
in the sense of Cauchy principal value, which may be calculated as in [TT] [T8] . We find 



(L-h-T), (29c) 

which are the modal equations that generalize those of Waterman and Truell, eqs. 
in [33]. The extinction theorem ( |2"7|) is unchanged. 



3 Matrix form of the modal equations 

In this Section we focus on reformulating either of the identical systems of equations (1261) 
with the objective of determining a single scalar equation for the wavenumbers. With no 
loss in generality we consider the first system, which in full is as follows: 
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The goal is a simplified equation for the wavenumbers £ and 

We begin by rewriting the doubly infinite set of equations in matrix form. The un- 
known amplitudes are represented by vectors a a , defined by a Q = (. . . , A°_ X) Aq, Af, . . .)*, 
a G {L,T}, so that the amplitude A^ is in the pth position of the infinitely long column 
vector. Introduce the constant vector e = 
trices Q° and T a/3 with elements 



, 1, 1, 1, 



Q 



in 

a _ _2_ 

ran 



N a (£) - 1 



^^rL, a,/3e{L,T}. 



Then eq. ( 1301) can be expressed 



*, and the infinite square ma- 

(31) 



a L - e(Q L + — ) (T"a L + T Ti a T ) = 0, 

Vl 

a T — e(Q T + — ) (T LT * L + T"a T ) = 0, 

Vt 



where the scalars e and y a are 



e = — i4n , 
2/« = e 2 -^, ae{L,T}. 



(32a) 
(32b) 



(33a) 
(33b) 



The pair of equations 
vectors and matrices 



with vectors 



and matrices 



are combined to form a single equation involving doubly infinite 

(34) 



{X - eQT - e^T - e^Tja = 

Vl Vt 



I 
I 



Q 



Q L 
Q T 



a L 

8lt 



rj-i.LT rj-iTT 



(35) 



(36) 



The system (134)) can be replaced by an equivalent equation of simpler form. In order to 
achieve this we multiply from the left by T 1//2 and rearrange the terms, 



{X - eQ - -f L fl - — f T f^}u = 0, 

Vl Vt 



where the doubly infinite vectors are 

f L = T^e L , f T = T 1 / 2 e T , 
and the sole matrix, apart from the identity, is 

Q = T^QT 1 / 2 . 



u 



T l/2. 



(37) 

(38) 
(39) 
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The wavenumbers must clearly satisfy the condition that detP = where P is the 
infinite matrix premultiplying u in eq. ( 1371) . This approach was adopted by Varadan et 
al. [30] , although the infinite system of equations they considered are not the same as the 
present system ( 1381) . While detP = is a sufficient condition for determining the coherent 
wavenumbers, it is not always necessary, particularly in the limit of small concentration, 
which is equivalent to the limit of small e. We now develop a necessary condition that is 
far simpler in form, and which, as we note below in §4.11 is both necessary and sufficient 
for small e. We first rewrite (|37|) in the equivalent form 



{X-e(X-eQ) L ( 



1 i f if L , 1 T 1 T 



Vl 



+ f ^)}u = 0. 



Vt 



(40) 



Taking the inner product of ( 1401) with i L and i l T yields a system of two equations for the 
scalar quantities f^u, a G {L,T}, 



1 - e 



Mll 
Vl 



■ M LT 
Vl 



The four matrix elements are given by 



- Mtl 
Vt 



Vt ■ 



'f*y 



(41) 



M. 



al3 



f* a (X-eQ) 1 f >9 , a,/3G{L,T}. 



(42) 



Taking the determinant of the 2x2 matrix yields the desired equation for £: 

(j/i - eM LL ) (y T - eM TT ) - e 2 M LT M TL = 0, (43) 

or more explicitly 

(44) 



(e -k\- eM LL ) (e - 4 - eM TT ) - e 2 M LT M TL = 0. 



Equation ( |44"|) is the fundamental equation for determining the coherent wavenumbers. 



4 Asymptotic solutions of the wavenumber equation 
4.1 Preliminary observations 

Before considering solutions to eq. (I44p we note that the matrix elements can be expressed 



M, 



= e^T 1/2 (X- eT 1/2 QT 1/2 ) _1 T 1/2 e /3 
= <[T- 1 -eQ]- 1 e /3 , «,/3e{L,T}. 



(45) 



This indicates that the elements can be calculated without evaluation of the square root 
matrix T 1 ^ 2 . The elements depend upon £, but the form of the wavenumber 
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equation (jUJ) suggests a natural asymptotic expansion in the parameter e. Thus, using 
(X - eQ)- 1 = X + eQ + e 2 Q 2 + . . ., eq. (02) implies 

M Q/3 = + eflQf^ + e 2 f* a Q% + . . . 

= e^Te^ + ee^TQTep + eV Q TQTQT e/3 + . . . , a, (3 e {L, T}. (46) 

Based on this expansion, the modal equation eq. (1441) has the following leading order 
expansion in n 



This indicates that even at low concentration the effective wave numbers which are so- 
lutions of eq. ( 143]) depend on a coupling between P and SV waves. It should also be 
noted that eq. ( 1471) is an algebraic equation of order 2 with regard to £ 2 . Consequently, 
as expected at low concentration, there are two solutions that propagate in the direction 
of the x— axis (the other ones propagate in the opposite direction). 

Satisfaction of the single equation ( T4*4"|) is clearly a necessary condition for the infinite 
system of homogeneous equations 034p . An alternative procedure to solving the latter 
for small e is to simultaneously seek both the wavenumber £ and the null vector a as 
asymptotic series in e. We do not present the details here, but it can be shown that the 
leading order solution yields a = e a +0(e) for a = L or T. This indicates that the effective 
wave solution can be characterized as either a quasi-L or a quasi-T wave. Furthermore, 
it may be shown that the solution for the wavenumber obtained by proceeding with the 
asymptotic series in e is equivalent to that of the simplified equation (|441) . Hence, the latter 
is both necessary and sufficient for developing the asymptotic wavenumber solution. 

The remainder of this Section considers asymptotic expansions of the solutions, valid 
in different limits: low concentration, low frequency, and high frequency, respectively. We 
begin with the small concentration expansion. 



4.2 Asymptotic expansion at low concentrations 

Rather than working with the wavenumber directly it is more convenient to expand the 
solutions of eq. ([431 about one of the two leading order solutions y a = 0, a G {L, T}. We 
choose to expand about the P-wave root yi = 0, although all of the results below apply 
to the other solution under the interchange L -h- T. At low concentration |e| 1 and we 
therefore assume a formal asymptotic expansion in e as follows: 

Vl = ey? + ?yf + e'yf + .... (48) 
Substituting y L from (14"8"|) into eq. (j43|) and noting that yr = Vl + k 2 L — k\, we obtain 

(y? - M LL + eyf + e*yf + . . . ) (k\ - k\ + e(y? - M TT ) + ?yf + . . . ) 

- eM LT M TL = 0. (49) 

The coefficients in the expansion (1481) follow by taking derivatives of eq. ( I49p with respect 
to e at e = 0. 



Conoir & Norris December 31, 2009 



13 



4.2.1 O(e ) 

The leading order term is found by setting e = 0, yielding 



Hence, 



v { l ] - m ll 



y^ = e t T LL e = f LL {Q). 



e=0 



0. 



LLi 



(50) 



(51) 



4.2.2 0(e 1 ) 

At the next order eq. (|49p implies 

(2) 



dM 



LL 



de 



e=0 



M LT M' 



TL\ £=0 



0. 



(52) 



The derivative can be found using the expansion ( )46|) . resulting in 

e*T Ti ee*T iT e 



T LT e 



p _ p 

T 1 \L 



(53) 



The quantities Q L and Q T in (1531) are evaluated at £ = ki,. Using the definition of Q T in 



eq. (I3"T|) we have 



yf = e'T LL Q L T Li e + 



111 ^trnTL-^rT ,, xri-iLT, 



-eY L N J (fe)T LJ e, 



2(^1 4; 

where the square matrices N°, Qq and for later use, Qg , are defined 

NU0 = N«_ n (0, Qo = Q a (K), Qo =Q a '(k a ), ae{L,T}. 
Expanding the function iV"(£) for small (£ — k a ) it is straightforward to derive 

Qomn D m _ ri {k Q ), Qomn D m _ n (k a ) , 

where 



(54) 



Df\k) 



111 

T 



{{kbf -p 2 )J p {kb)H^\kb) + (kb) 2 J' p (kb)H^\kb) 
p 2 {kb) 2 



l + i2^(kb) 2 J p (kb)H ( p l \kb) 



(55) 

(56) 

(57a) 
(57b) 
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4.2.3 0(e 2 ) 

Taking the second derivative of ( 149]) and setting e = yields 



e=0 



2[yf _ dA^ 



e=0 



] [y™ - M TT | e=0 ] - ^-MltMtl 



= 0. (5£ 

e=0 



The derivatives are again evaluated using the expansion ff46|) . The second derivative re- 
quires some care, since the terms involving Q on the right hand side of f l46]) are themselves 
dependent upon e. This is the same as for the acoustic case; in fact, in the acoustic case 
only the second derivative terms is needed, all the remaining terms in f )58|) are unique to 
the elastic problem. Thus, using results from the acoustic case we have 

(f TT (0) - f LL (0)) 

yf = eiTQTQTe L + / LL (0)e i L TQ / Te L + ^ _ fc 2 )2 ^ e' L Ter e£,Te £ 

+ 2 (kl-k 2 T ) [ (e ^ TeT) e TTQTe L + (e^Te L ) e^TQTe T ] , (59) 

where the doubly infinite square matrices Q and Q' are evaluated at £ = k^. Higher order 
terms in the expansion fT4*8l) can be determined using the same procedure. 



4.2.4 Wavenumber expansion to 0(nfy 

Combining the above expressions for y^' and the quasi-longitudinal wavenumber 
expansion up to second order in the concentration is 



C = ki + dfn + <%r4 + . . . , (60) 



where 



d!i = -4i/^(0), (61a) 
4 = -SE - E ^M^lf T . (61b) 

The third order term follows in a straightforward manner from y L but is too long to 
warrant including here. The expansion for the quasi-transverse wavenumber follows by 
the interchange L -h- T '. 

4.3 Long wavelength limit 

The long wavelength regime corresponds to small k^^b. As k a b — > we have 

W^UhY- ae[L ' n (62) 
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In the same limit, it follows from its definition in eqs. ( 13~TT) and ( |55i) that the elements of 
the matrix Qq become 

Qomn = \\ m ~n\, (63) 
The long wavelength limit of the low concentration expansion ( )60l) becomes, up to O(np), 

e = k 2 L + 5^n + + . . • , (64) 

where 

if = -4i/"(0), (65a) 



Note that the expansion for the wavenumber is independent of b in this limit. Again, the 
quasi-transverse wavenumber expansion follows from the interchange L -h- T. 



4.4 High /cl,t& limit and the Waterman & Truell formula for 
elastic media 

As mentioned in the Introduction, the method we have developed is also valid at high 
frequency. In the following it is assumed that k^^b — > oo. This means either the frequency 
increases, i.e. ki^TO, — > oo with b = 0(1), or the separation length increases, b — > oo with 
^l,t<2 = 0(1). The high frequency limit k^^a oo corresponds to wavelengths shorter 
than the radius a of cylinders. This situation has been studied in [6] which deals with 
ultrasonic characterization of thermal damage in concrete. However, since b can be related 
to the concentration c (0 < c < 1) of cylinders [3lj by 



n na 2 a 2 



riQirb 2 b 2 



(66) 



it follows that the limit b — > oo with fc^a = 0(1) corresponds to a dilute medium 
(c — > 0). This is one reason why it is of interest to compare the high k^^b limit to the 
Waterman & Truell approach. 

In the limit as k^^b — > oo the leading order asymptotic approximation of iV^(£) is [1] 

iV p a (0 * ^^/y [(-l) p (*a -Oe* 6 +»(*a + 0e"* 6 ] e^ b , « = {L,T}, (67) 

It follows that the coupled systems of infinite equations ( 1301) become in this limit 

A^ + P LT + (-l) n Q L T = 0, (68a) 
(L-H-T), (68b) 



CONOIR & NORRIS 



December 31, 2009 



16 



with 



LT 



p 

Qlt = Ql ^2(-l) p (Tp L Ap + Tp L Ap), 
p 

(L^T), 



and 



The structure of eqs 




Q< 




( )68l) implies the identities 



A_ n — A^' T and A n ' +2 



2n I k a e l ( ka+ ^ b 



a 



A 



L.T 



(69a) 
(69b) 
(69c) 

{L,T}. (70) 
(71) 



Consequently, the problem reduces from calculating an infinite set of unknowns to one with 

L T 

eight unknowns: Plt, Ptl, Qlt, Qtl and A \ satisfying a system of eight homogeneous 

linear equations. Note that although Plt, ■ ■ ■ ,Qtl can be expressed in terms of A \ , 
the calculations are simpler with Plt, - ■ ■ ,Qtl considered as unknowns. The first four 
equations are obtained by setting n = and n = 1 in eqs. (|68|) . Then, we perform an 
iteration on eqs. (1681) using the identities (17TT) . with the result 



A L n - [f LL (0)P L + (-irf LL (<K)Q L ]P LT - [f LL (-K)P L + (-lTf LL (0)Q L }Q L T 
- [f TL (0)P L + {-l) n f TL {7t)Q L ]P TL - [f TL (n)P L + (-lTf TL (0)Q L ]Q TL 



(L^T). 



0, 

(72a) 
(72b) 



Four equations are obtained from eqs. ( 1721) by considering the two possibilities for (— l) n , 
corresponding to n = and n = 1. The second set of four equations follow from eqs. 
( 1691) combined with the identities (ITTi) . Eliminating Aq^ from the eight equations results 
in the following four equations for the four unknowns Plt, Ptl, Qlt, Qtl, 



[i + f LL (0)P L + f LL (*)Q L ]PLT + [i + f LL (*)P L + f LL (0)Q L ]Q 



LT 



[f TL (0)P L + f TL (n)Q L }P T L + [f TL WP L + / TL (0)Ql] Qtl = 0, (73a) 



[1 + f LL (0)P L - f LL (n)Q L ]P L T - [i - f LL ^)P L + f LL (o)Q L ]Q 



LT 



+ [f TL (0)P L - f TL (n)Q L }PTL + [f L {iz)P L - f TL (0)Q L ]Q T L = 0, 

(L^T). 



(73b) 
(73c) 

The homogeneous linear system of equations f!73|) has nontrivial solutions if the associated 
determinant vanishes. Thus, the modal equation at high frequency is 



det 



1 + f LL (0)P L 

f LL ^)Q L 

f LT (0)P T 
f LT (n)Q q 



*LLi 



>T 



f LL ( 

1 + / LL (0)Ql 
f LT WP T 

f LT (o)Q T 



f IL (o)P L 
f TL (*)Q L 

i + f TT (0)P T 
f TT ('k)Qt 



f 1L {*)p L 

f TL (o)Q L 

f TT (7T)P T 

i + f TT (o)Q 



T 



(74) 
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The mode converted forward scattering and back-scattering amplitudes, f LT,TL {0) 
and f LT ' TL (ii), respectively, are identically zero if the fundamental scatterer has sufficient 
geometrical symmetry. This is the case for circular cylinders, and occurs generally for 
cylinders with reflection symmetry about the x— axis. When f LT ' TL (0) = f LT ' TL {ji) = 0, 
instead of the determinant of eq. (I74p . the condition for satisfaction of the four equations 
(173"j) becomes two simpler equations: 



[1 + f LL (0)P L ] [l + f LL (0)Q L ] - [f LL (n)] 2 P L Q L = 0, 
[1 + f TT (0)P T ] [l + f TT (0)Q T ] - [f TT (7r)] 2 P T Q T = 0. 



(75a) 
(75b) 



These provide uncoupled modal equations for the P and SV waves, in contrast to what 
happens with the Fikioris & Waterman approach. 

Starting from eqs. (I29j) and using the same process as before, we can easily verify that 
the Waterman & Truell's modal equation for elastic media is once again given by eqs. 
(174)) where now P^t and Ql,t are defined by 



P n 



2nc 



Q, 



2n 1 



a 



{L,T}. 



(76) 



So, the modal equation in the limit k^^h — > oo is very similar but fundamentally different 
from that of Waterman & Truell. For circular cylinders, eqs. (1751) combined with ( 1761) 
are nothing else than the Waterman & Truell formula [33] for acoustic waves applied to 
P and SV waves separately. 



5 Reflection coefficients 

5.1 Reflection coefficients for an incident P wave 

We consider a single type of wave incident on the half-space x > 0, specifically a P wave, 
so that Al — 1 and At = 0. The fields {ip L {r)) and (tp T (f)) in the domain {x < 0} then 
correspond to the P and SV reflected waves respectively. Consequently, the reflection 
coefficients R LL and R LT can be defined as follows 

(ifj a (r)} = R La e~ ikaX , (x<0), a = {L,T}. (77) 

Inserting eqs. ( ITU]) and (fTTj) into eqs. ([7]), and taking into account eqs. (fTTj) . we get 



<^(r)> = n J2 i~ P (T p LL A L p + TjMJ) jft) + (T^B L p + T^fiJ) J^) 
p 



with (fji = (xji^jx)) 



J?{0=noe* x f df^Kr^e^, ( E {£,£'}, a E {L, T}. 

J Xjl>—X 



(78a) 

(78b) 
(78c) 

(79) 
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The integrals in eqs. ( |79l) are calculated in (TTJ [18], with the result that 
and consequently 



o- / r pLL AL i rpTL AT t>LL tdL _j_ t-iTL dT 

*» - £ B- < ' * + ) , (SXa) 



o- /rnLT AL , rpTT AT rpLT r>L _i_ rpTT dT \ 

^ LT = 2mo / ip p p A> + fg ^ ) . (8 i b) 



V 

As expected, R LL and i? iT depend on the infinite sets of coefficients A^ T and B^ T . 
These coefficients follow from eqs. ([27)1 . which we rewrite in matrix format: 

where the notation is the same as in Section [31 with the additional quantities b* = 
(t4b*0, b Q = (..., B Q , Sf ,...)*, and ^ = £' 2 {L,T}. The infinite 

vectors a and b are null vectors of (|34j) corresponding to the quasi-P and -SV roots of eq. 
(j44j) . respectively. To be specific, {£, a} solve (|34j) and {£',b} satisfy 

{X-eQT-e^^T-e^^T}b = 0, (83) 
Vl Vt 

which together determine the null vectors up to scalar multiples. Finally, eqs. ( 18 2 p fix the 
amplitudes of the null vectors in terms of the incident wave amplitudes (Al = 1, A? = 0). 

Once the coefficient vectors a and b are determined, the reflection coefficients follow 
from (IHTj) . which can be recast in matrix form as 

R LL = e (^-0 e y Ta + e ( k *<-p e tj Th (84a) 

2k L y L 2k L y' L 1 ; 

R LT = £ %^e^JTa + e^^e^JTb, (84b) 
2k T y T 2k T y' T 1 ; 

where J = diag(. . . , —1, 1, —1, 1, . . .) with J o — 1- The reflection coefficients R TT and 
R TL are obtained in the same same way with Al = and At = 1 in eqs. 



5.2 Asymptotic approximation 

As in Section [31 we consider the low concentration asymptotic approximation. Starting 
with the ansatz 

a = a (°> + ea« + . . . , (85a) 
b = b (0) + eb (1) + . . . , (85b) 
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and using the result from Section H] that yi = ee^Te£+0(e 2 ) we find that, to leading order 
a (°) = A L e L . Similarly, y' T = ee T Te T +0(e 2 ) and hence b (0) = A T e T = 0. Equations (E 
then become 



R LL = + 0(e 2 ), (86a) 



* 1T — 2TOS + < £2 >' <86b) 

The case of SV wave incidence (Al = 0, At = 1) can be treated in the same way. 
Combined with the previous results, we find that the leading order approximations to the 
reflection coefficients are: 

^=^tS +0( "» ) ' a '^^- < 87 > 

This result shows that reflected waves do not depend on the effective wave numbers at 
this order, and helps explain why it can be difficult to detect reflected coherent waves. 
Experiments usually deal with transmitted coherent waves [8]. 



6 Generalization of the Linton & Martin formula 

The aim of this section is to generalize the acoustics formula of Linton & Martin [T8j 
to elasticity. The advantage of the Linton & Martin formula is that it expresses the 
wavenumber at low concentration in terms of the far-field scattering function only, rather 
than the T-matrix elements. In the present context of elastodynamics, this requires that 
we express the two series in (165bl) as integrals of the far-field scattering functions. The 
analogous acoustic problem involves only the first series in (165b|) , which was calculated in 
[IE] as 

- 4 E I- "I W = ^ jfd S «*$±U" m >. (88) 

L m,n L ,/u 

This integral is, however, converg ent only if f LL (0) = 0. That is the case for certain 

scatterers, including circular cylinders, where T^h = T^ L =3- f LL (0) = 0. But the 
formula is not correct for non circular cylinders and must therefore be modified. The 
resolution of this problem is found in the following identity, which follows from eq. ( JT3|) . 

f LL (0)f LL (-0) = J2 T n LT m L cos(n - m)0. (89) 

m,n 

Using this we derive an identity similar to eq. (83) of [TS] that does not rely upon the 
assumption = T^ L . Thus, instead of (}8"8"]) we have the more general identity 

- ^ E I™ - "l T " T n L = ~j^2 f^^ ^t( d -)^[f LL (9)f LL (-9)}. (90) 

L m,n L ^° 
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Regarding the second series in ( I65bl) and its dual for the quasi-SV wavenumber, we 
consider 

^ £ ^ T ™ T n L = \ £ ( W + Ola) 

m,n m,n 

^ £ (|) Im " n ' T ^ T n T = \ £ + T ^ T n T ) ^ (91b) 

m,n m,n 

where k = hx/k^ > 1. Although it might appear that the series for St is divergent, 
it is not. Based on the original result of Muller [23] for the Helmholtz equation, later 
extended to elasticity [7], it follows that the far-field scattering functions are necessarily 
entire functions of the angular argument 9 considered as a complex variable. In the present 
context this implies that the Fourier coefficients of the far-field scattering functions decay 
in such a manner that St series is convergent for any finite k. This analytic property 
is apparent by considering, for example, the far-field scattering functions for a circular 
object. 

Our main result is that the two sums can be simplified as follows: 
Lemma 1 

o _ ^ - 1) r d9G LT (0) (q2) 
L 2ir L 1-2kcos6 + k 2 ' { ' 



and 



S L + S T = Gurii log k) 



lim ^ - / d6G LT (9)1 - '— \k n \, 93 

tv^ooItt / LM J \ 1-2kcos# + k 2 



G lt (0) = f" {6)f 1L {-6) + f lh {6)f L1 (-9) = G LT {-6). (94) 



where 



Note that the limit in fl93l) does not commute with the integral. This is a consequence of 
the properties of the far-field scattering functions mentioned before. 

Proof: Starting from the definition in eq. ( 19~4"|) and using eq. f JT3|) . we have 



Glt{B) = £[W + T T m L T^\ cos(m - n)6. (95) 

m,n 

Consider the function 

g (0) = y K - |n| cos^ «rH = i/ d9g(9) cos n9. (96) 

^ Wo 

n 

Hence, eq. fl95|) implies that 

S L =^j\0g{e)G LT {9). (97) 
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Performing the sum ( l96l) i yields 



from which the result (l92j) follows. It remains to prove (|93l) . 
The first identity in (193]) follows directly from ( |95|) and 

Glt(* log «) = \ Yy^l L + Tl L T^ T ] («-» + «»-») . (99) 

m,n 

Regarding the second identity in f l9"B"j) . note that replacing k with in eq. fl9"oT) i gives 
a divergent series. However, the following finite series is perfectly well defined for any 
positive integer N, 

h N {9) = V K H cosn0 K lnl = -[ d6h N (6)cosn6, \n\ < N, (100) 

\n\<N 71 J ° 

and therefore 

£ ^ T ^ + ^^ (101) 



m. n 

|m — n| < iV 



Using 1 + k + k 2 + . . . + = (1 - k n )/{1 - k) gives 

. /.x ~ N f kcos(N -1)9 -cos N9\ 

which implies the second identity in ( |93i) . The integral formula ( 1931) provides a formal 
procedure to evaluate the far-field scattering functions at an argument that is imaginary, 
such as ±i log k in this case. In practice, it is advisable to use the Fourier series, or 
perhaps a truncated version of the same if there is any numerical noise present. This is 
equivalent to a singular value decomposition (SVD) of the far-field scattering operators, 
which is always necessary in any inverse scattering algorithm. 

Combining the simplified expressions for Sl and St with the expansion f|6%|) . we arrive 
at our main results: 

Theorem 1 The long wavelength limit of the low concentration expansions for the quasi- 
P and -SV wavenumbers are, to 0{nfy, 

2tt 



e = kl - 4m /"(0) + ^J'dB cot( 9 -)±[f LL (9)f LL 



K }\ R f LT (o)f TL (-o) + f TL (Q)f LT (-o) rin ^ 

7t Jo k T -2k L k T cos9 + k 2 L ' 1 ' 
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and 

e = 4 -4m o/ ^(0) + ^ jfd* cot( 6 -)-^[f TT (9)f TT (-9)} 

+ 8nf r dg / LT w/ TL (-^) + / TL w/ Lr (-^) 

7r 7 fcf, — 2/c L fc T cos 9 + k\ 




(104) 

These formulae generalize the identity found by Linton & Martin [TS] for acoustic 
(scalar) waves. The acoustic equation can be derived from f 1 1 U 3 j) by eliminating the SV 
wave contributions associated with the index T, giving 

e = kl-4 ino f LL (0) + ^ jTd* cot( 6 -)f- e [f LL (9)f LL (-9)]. (105) 
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